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ABSTRACT 

The  Aerospace  Nuclear  Safety  Program  at  Sandia  Corporation  requires  an  independent 
assessment  of  the  safety  of  all  nuclear  power  supplies  used  in  aerospace  under  all  conditions 
from  factory  shipment  through  re-entry  from  space.  The  effects  of  impact  against  any  point 
on  the  earth  by  a  nuclear  power  supply  or  its  components  must  be  evaluated  as  part  of  this 
assessment.  This  report:  (1)  evaluates  the  fundamental  mechanisms  involved  in  soil  impacts 
by  cylindrical  projectiles  with  velocities  of  500-700  feet  per  second;  (2)  defines  the  problem 
of  specifying  impact  targets  which  closely  simulate  various  portions  of  the  earth's  surface; 
and  (3)  describes  an  analytical  and  experimental  program  intended  to  obtain  the  information 
needed  to  properly  specify  a  simulated  earth  target  for  use  in  Aerospace  Nuclear  Safety 
impact  tests. 
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PRELIMINARY  EVALUATION  OF  EARTH  TARGETS  FOR  USE  IN 
IMPACT  EFFECTS  STUDIES 


I.  Introduction 


The  Aerospace  Nuclear  Safety  program  at  Sandia  Corporation  requires  an  independent  assessment 
of  the  safety  of  all  nuclear  power  supplies  which  are  to  be  used  in  aerospace  applications.  This  safety 
assessment  must  cover  all  credible  hazards  which  might  occur  under  any  conditions  that  could  exist  from 
the  time  of  factory  shipment  through  re-entry  from  space.  One  major  condition  which  requires  analytical 
and  experimental  investigation  is  that  of  power  supplies  and/or  their  components  impacting  any  point  on  the 
surface  of  the  earth  at  varying  velocities. 

In  order  to  provide  data  and  information  of  value  for  assessing  the  safety  of  nuclear  aerospace  power 
supplies,  the  ground  tests  must  be  cogent,  engineering-type  tests  which  are  so  well  planned  and  performed 
that  test  conditions  can  be  closely  controlled  to  simulate  actual  conditions.  In  the  case  of  earth  impact 
tests,  conditions  of  velocity,  temperature,  attitude,  and  structural  strength  of  the  test  item  at  impact  must 
be  closely  controlled.  These  requirements  obviously  preclude  actually  dropping  a  test  item  from  an  air¬ 
craft  against  the  natural  earth  because  it  is  impossible,  practically,  to  sufficiently  control  the  conditions 
of  test  item  attitude,  temperature,  and  velocity  during  a  free-fall  drop  test  to  produce  the  required  cogent 
data.  The  only  recourse  remaining  is  to  conduct  earth  impact  tests  at  a  test  laboratory  location  where  the 
test  conditions  can  be  closely  controlled  and  to  use  a  test  target  that  closely  simulates  several  representa¬ 
tive  areas  of  the  earth's  surface. 

A  survey  of  impact  targets  used  in  previous  tests  to  simulate  the  earth's  surface  reveals  that  these 
targets  have  not  closely  simulated  the  earth's  surface.  Recognizing  this  situation,  Sandia  Corporation 
has  launched  a  program  to  evaluate  the  requirements  of  an  earth  target  and  to  develop  a  specification  for 
a  series  of  such  targets.  This  report  discusses  this  program. 


II.  Statement  of  Program  Objective 


The  program  described  in  this  report  will  evaluate  the  requirements  of  an  earth  impact  target  and 
develop  a  specification  for  a  series  of  impact  targets  which  will  adequately  simulate  the  variety  of  earth 
surfaces  that  might  be  impacted  by  a  nuclear  power  supply  intended  for  aerospace  applications. 

The  attainment  of  this  objective  implies  a  program  consisting  of  the  following  phases: 

1.  A  general,  analytical  study  of  fundamental  mechanisms  involved  in  the  penetration  of  earth 
media  (ranging  from  solid  rock  to  water)  by  cylindrical  objects  traveling  at  500-700  feet 
per  second. 

2.  Studies  of  model  projectile  penetrations  into  various  media  to  empirically  relate  these  laboratory 
results  to  the  results  of  the  analytical  study. 

3.  Full-scale  studies  of  vertical  projectile  penetrations  into  various  simulated  earth  media  to 
empirically  relate  the  results  observed  to  the  results  of  previous  model  studies. 

4.  A  preliminary,  full-scale  study  of  horizontal  projectile  penetrations  into  various  simulated 
earth  targets  to  empirically  relate  results  observed  to  those  of  previous  full-scale,  vertical 
penetration  studies. 
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5.  Preparing  a  tentative  specification  for  a  simulated  earth  target  which  will  be  used  for 
horizontal  rocket  sled  tests. 

6.  An  experimental,  full-scale  study  of  projectiles  penetrating  horizontally  into  various 
simulated  earth  targets  which  are  fabricated  to  the  previously  prepared  tentative 
specification. 

7.  Preparing  a  final  specification  which  defines  earth  targets  for  horizontal  rocket  sled  tests 
which  closely  simulate  actual  earth  impacts. 

Various  corollary  studies  (such  as  model  similitude,  water  impact  mechanisms,  earth  viscosity 
variations,  and  compressibility  and  stress-strain-time  relationships  for  various  soils  under  varying 
rates  of  strain)  will  also  be  performed  to  help  achieve  the  objectives  of  this  program. 


III.  The  Earth's  Surface  as  a  Target  Material 


It  is  possible  for  an  aerospace  nuclear  power  supply  or  some  of  its  components  to  impact  upon  any 
part  of  the  earth's  surface  during  certain  conditions  of  launching  and  of  flight  duration.  To  completely 
assess  the  safety  of  these  aerospace  units,  the  effects  of  such  impacts  upon  the  earth's  surface  must  be 
evaluated.  Ground  impact  tests  to  evaluate  these  effects  require  targets  which  simulate  representative 
areas  of  the  earth's  surface. 

The  surface  of  the  earth  is  far  from  being  a  homogeneous  and  isotropic  medium.  It  ranges  from 
very  hard  exposed  rock  through  many  varieties  of  soil  cover  to  water.  In  a  very  few  feet,  the  properties 
of  the  solid  portions  of  the  earth's  surface  can  vary  greatly  in  both  horizontal  and  vertical  directions. 
Although  soil  maps  and  geological  maps  have  been  made  of  many  parts  of  the  surface  of  the  earth,  much 
of  the  total  land  area  has  not  been  accurately  mapped.  The  U.  S.  Corps  of  Engineers  and  other  agencies 
arc  currently  mapping  uncharted  surface  areas  all  over  the  world.  When  and  if  this  mapping  is  completed, 
it  may  be  possible  to  statistically  predict  the  probability  of  an  aerospace  power  unit  impacting  a  particular 
typo  of  surface.  However,  in  order  to  evaluate  the  effect  of  a  power  plant  impacting  a  given  type  of  surface 
it  will  be  necessary  to  perform  impact  tests  against  a  limited  number  of  simulated  earth  targets  which  have 
been  arbitrarily  selected  as  typical  of  the  various  kinds  of  earth  surfaces. 

Although  it  is  possible  to  classify  the  surface  of  the  earth  as  rock,  soil,  or  water,  in  nature  these 
materials  arc  blended  to  form  a  varying  medium.  Even  uniform  media  exhibit  variations  in  strength. 

For  example,  the  harder  igneous  rocks,  such  as  basalt  and  felsite,  have  the  highest  potential  shear 
strength,  but  even  these  rocks  have  joints  and  cracks  which  weaken  the  material  non-uniformly. 

Soil  is  composed  of  water,  air,  and  mineral  particles.  Its  strength  and  in  situ  modulus  of  deforma¬ 
tion  depend  on:  (1)  the  type,  shape,  and  size  of  the  mineral  particles;  (2)  the  structural  arrangement  of 
the  mineral  particles;  (3)  the  moisture  content;  (4)  the  density;  (5)  the  original  confining  pressure;  (6)  the 
time  rate  of  loading;  and  (7)  the  stresses  caused  by  the  loading  element.  The  ultimate  strength  of  a  soil 
also  depends  upon  the  size  of  the  loading  clement  and  is  highly  time  dependent. 

Soils  range  from  sands  and  gravels  which  have  low  void  ratios  (the  ratio  of  the  volume  of  the  voids 
in  a  soil  to  the  volume  of  the  solids)  to  very  fine,  sensitive  clays  with  high  void  ratios  whose  moisture 
content  may  vary  as  much  as  several  hundred  percent.  In  most  soils,  the  shear  strength  increases  as  the 
void  ratio  and  moisture  content  are  reduced;  however,  this  relationship  is  not  always  true. 


No  such  thing  as  a  "typical  soil"  exists  in  nature.  Therefore,  one  must  arbitrarily  select  a  limited 
number  of  soil  types  to  be  simulated  so  that  the  effects  of  impact  upon  any  soil  can  be  reasonably  predicted. 
Unfortunately,  the  proper  soil  types  for  simulation  cannot  yet  be  selected  without  more  information. 
However,  the  acquisition  of  this  information  will  be  one  of  the  results  of  this  entire  program.  To  limit 
soil  variations  for  better  understanding,  the  earth's  surface  might  be  classified  into  several  arbitrary 
groupings  of  varying  shear  strength. 
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Rock 


Sand 


Silt 


Clays 


1.  Massive  igneous  rock  -  basalt,  diabase,  felsite,  etc. 

2.  Stratified  limestone 

3.  Coarse  gravel  and  fragmented  material 

4.  Dense  sand  with  low  moisture  content 

5.  Loose  blow  sands 

6.  Saturated  dense  sands 

7.  Dry  cemented  loose  silts  or  loess 

8.  Saturated  silts 

9.  Desiccated  clays 

10.  Overconsolidated  saturated  clay 

11.  Normally  consolidated  saturated  clays 

12.  Water 


NOTE:  These  groups  are  not  necessarily  arranged  in  order  of  descending 
resistance  to  penetration  or  shear  strength,  but  are  arranged  in 
order  of  descending  size  of  the  particles  making  up  the  media. 


In  each  of  the  above  groups,  the  lateral  and  vertical  pressures  increase  with  depth.  For  earth 
pressure  at  rest,  there  are  no  shear  stresses  in  the  horizontal  and  vertical  planes.  By  definition,  the 
value  of  the  lateral  earth  pressure  existing  in  the  undeformed  mass  of  natural  soil  with  a  level  surface 
is  called  the  "lateral  earth  pressure  at  rest.  " 


As  can  be  seen,  the  earth's  surface  is  a  complex  target  presenting  media  of  varying  strength  to  the 
impact  of  an  incoming  projectile. 


This  report  discusses  the  mechanism  of  penetration  in  the  various  media  outlined  above,  in  order 
to  develop  criteria  for  selecting  a  limit  number  of  "typical"  soils  to  study  the  impact  of  nuclear  power 
supplies  in  free  fall  on  the  earth's  surface. 


IV.  Simulated  Earth  Targets  Presently  Being  Used 


A.  Description  of  Present  Targets 


For  many  years,  simulated  earth  targets  have  been  used  in  horizontal  rocket  sled  impact  tests  at 
many  laboratory  locations.  In  fact,  several  impact  tests  of  aerospace  nuclear  power  supplies  have  been 
performed  using  simulated  earth  targets.  As  explained  earlier,  the  use  of  horizontal  rocket  sled  tracks 
is  highly  desirable  for  these  tests  because  a  variety  of  test  conditions  can  easily  be  controlled  closely  to 
insure  meaningful  results. 

Earth  targets  at  Sandia  Corporation  consist  of  compacted  sandy  silt  with  a  small,  undetermined 
amount  of  colloidal  clay.  The  material  is  excavated  from  its  natural  site  near  the  Sandia  sled  track. 

After  being  moistened,  it  is  compacted  in  approximately  6-inch-thick  layers  by  a  hand-held  air-powered 
compactor  into  an  8-foot  x  8-foot  x  8-foot  wooden  form  at  approximately  95  percent  of  Standard  A.  A.S.H.O. 
Density.  The  unit  weight  of  the  excavated  material  and  the  moisture  content  of  the  compact  material  are 
unknown. 

The  dimensions  of  the  target  are  selected  so  that  the  test  item,  after  impact,  is  completely  embedded 
and  has  come  to  rest  before  the  shock  wave  it  initiates  reaches  the  outer  surface  of  the  target  and  disinte¬ 
grates  it.  High-speed  motion  picture  films  of  past  tests  reveal  that  the  test  item  does  come  to  rest  before 
the  target  is  disrupted. 


As  presently  used,  the  earth  target  is  free-standing  above  the  surface  of  the  ground.  The  sides 
arc  not  confined  or  restrained  in  any  manner.  No  attempt  is  made  to  reproduce  the  normally  in  situ 
stressed  condition  of  the  soil  when  it  is  remolded  into  the  target. 

B.  Possible  Inadequacies  of  Present  Targets 

The  ultimate  shear  strength  of  a  soil  depends  on:  (1)  the  constitutive  properties,  (2)  the  effective 
mean  pressure,  and  (3)  the  rate  of  strain.  The  constitutive  properties  include:  (a)  particle  type,  size, 
and  shape;  (b)  soil  structure  and  degree  of  cementation;  (c)  void  ratio;  and  (d)  moisture  content.  The 
effective  mean  pressure  depends  on:  (a)  the  boundary  conditions,  (b)  initial  porewater  pressure,  (c) 
stresses  caused  by  the  penetrating  item,  (d)  the  rate  of  drainage  and  (e)  stress  history. 

The  total  lateral  pressure  of  the  earth  at  rest  is  very  difficult  to  measure,  but  it  is  approximately 
equal  to  0.4  to  0.  7  times  the  vertical  pressure  which  is  the  total  weight  of  the  material  above  the  point 
being  considered.  Effective  pressure  is  equal  to  the  total  pressure  minus  the  porewater  pressure.  The 
porewater  pressure  at  a  point  acts  equally  in  all  directions.  In  the  case  of  normally  consolidated  plastic 
clays,  the  lateral  earth  pressure  might  be  very  nearly  equal  to  the  vertical  pressure  because  of  creep 
and  stress  relaxation  with  time. 

If  the  surface  being  considered  is  at  rest,  there  are  no  shear  stresses  on  either  the  vertical  or 
horizontal  planes  and  the  vertical  and  lateral  pressure  increases  with  depth. 

Present  earth  targets  have  no  external  confining  pressure  increasing  with  depth,  but  are  held  in  place 
solely  by  internal  forces  or  cohesion.  It  is  the  cohesive  strength  (which  is  a  function  of:  (1)  moisture 
content,  (2)  particle  type  and  orientation,  (3)  void  ratio,  and  (4)  stress  history)  that  allows  the  present 
target  to  stand  vertically.  If  the  target  were  constructed  to  some  greater  height,  it  would  fail  in  shear 
because  of  its  own  weight  once  the  forms  were  removed.  If  the  target  material  had  only  frictional  strength, 
such  as  a  dry  or  saturated  sand,  it  would  fail  in  shear  once  the  forms  were  removed  regardless  of  the 
height.  In  that  case,  the  angle  of  inclination  of  the  faces  of  the  cube  would  approach  the  angle  of  internal 
friction  of  the  material.  After  the  forms  are  removed,  shear  resistance  is  mobilized  along  the  base  of  the 
cube.  The  weight  of  the  material  causes  the  sides  to  bulge  and  the  top  of  the  cube  to  subside.  This  process 
continues  until  enough  cohesive  and  frictional  resistance  is  mobilized  within  the  mass  to  keep  the  cube  in 
equilibrium.  The  cube  is  strained  throughout,  but  the  strains  are  most  pronounced  in  areas  behind  the 
faces  of  the  cube.  Tin's  is  usually  referred  to  as  the  slope  stability  problem.  The  soil  target  cube  is  in 
this  strained  state,  with  much  of  its  shearing  resistance  already  used  up  before  the  projectile  begins  its 
penetration. 

Therefore,  it  can  be  concluded  that  the  boundary  conditions  of  the  present  targets  are  in  no  way  the 
same  as  those  existing  in  the  same  material  located  in  situ  in  nature.  Likewise,  the  modulus  of  deforma¬ 
tion  and  other  fundamental  soil  parameters  are  not  the  same,  because  their  values  depend  upon  the  external 
boundary  conditions.  One  could  not  expect  the  results  from  impact  tests  performed  with  these  targets 
to  correspond  to  the  results  achieved  from  impacting  a  test  item  on  identical  material  in  its  natural 
location,  unless  the  boundary  conditions  in  the  target  are  made  to  equal  those  of  the  earth  at  rest. 

In  addition,  if  the  materials  used  in  the  target  were  compacted  to  exactly  the  same  density  and 
moisture  content  as  the  natural  soil  in  its  unexcavated  location  (a  situation  which  is  far  from  true)  these 
materials  still  would  not  have  the  same  shear  strength,  even  though  the  boundary  conditions  had  been 
identical.  This  is  because  the  excavated  material  has  been  extensively  remolded  and  its  structure  greatly 
changed.  It  is,  in  effect,  a  completely  different  material  after  disturbance  from  that  material  which  existed 
before  it  was  excavated  from  its  natural  location. 


V.  Hypothesis  of  a  Penetrating  Mechanism 


A.  General 

This  section  of  the  report  presents  a  hypothesis  of  a  possible  mechanism  for  a  cylindrical  projectile 
penetrating  into  an  earth  medium.  For  this  hypothesis,  the  example  of  the  impact  of  a  SNAP  10A  core 
vessel  as  a  projectile  is  used. 
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At  the  moment  of  impact  the  projectile  is  assumed  to  be  traveling  at  terminal  velocity  in  the  air 
with  no  horizontal  velocity  component.  The  projectile  weighs  150  pounds  and  the  impact  velocity  is  550 
feet  per  sec;  the  diameter  is  10  inches  and  the  length  15  inches.  The  total  energy  (E^)  is  equal  to  the 
sum  of  the  kinetic  energy  (EK)  plus  the  potential  energy  (Ep).  The  total  energy  that  must  be  dissipated 
before  the  projectile  comes  to  rest  is: 


=  ek  +  ep 


MV2 

2 


+  Wd, 


where 

M  is  the  mass  of  the  projectile  =  — -^0  lb.  force -  _  g6  slugs 

32.21  ft/sec/sec 

W  is  the  weight  of  the  projectile  =  150  lb.  force 
V  is  the  velocity  at  impact  =  550  ft/sec 
d  is  the  depth  of  penetration  in  feet. 

Thus 

Et  =  (4.66  s lugs M5 5 O' /sec)2  +  (15Q  lbs_  >  (d  ft) 

=  704,  362  ft.  lbs.  +  150  (d)  ft.  lbs. 

B.  Non- Penetrating  Projectile 


In  the  case  where  the  projectile  does  not  penetrate  the  surface,  rebounds,  or  is  destroyed  upon 
impact,  the  energy  will  be  dissipated  in  several  ways,  by: 

1.  Sudden  generation  of  heat  such  as  a  flash  of  light  and  sparks.  This  energy  loss  may  be 
insignificant  but  triboluminescence  and  frictional  heat  might  be  studied  to  evaluate  the 
magnitude  of  this  energy  dissipation. 

2.  Generation  of  sound  wave  caused  by  the  impacting  projectile.  This  energy  loss  is  probably 
insignificant  also,  but  some  study  can  be  made  of  its  magnitude. 

3.  Rise  in  temperature  of  the  projectile  and  the  soil.  There  will  be  a  loss  of  energy  because 
of  the  compression  of  the  projectile  and  the  compression  of  the  soil  surface.  Some  guess 
can  be  made  about  the  pressures  involved,  and  the  magnitude  of  this  energy  can  be 
estimated. 

The  heat  loss  is  probably  not  recoverable,  but  some  of  the  energy  used  to  compress  the 
soil  and  projectile  will  be  recovered  if  the  projectile  does  not  penetrate,  and  this  energy 
will  then  accelerate  the  projectile  back  into  the  air.  This  same  energy  may  also  deform, 
rupture,  or  fracture  the  projectile. 

4.  Compression  of  the  projectile  and  soil  will  cause  some  rise  in  temperature  of  both  soil  and 
projectile  and  is  not  recoverable.  Some  of  the  energy  used  to  compress  soil  and  projectile 
will  produce  the  force  to  rebound  the  projectile  or  its  parts  into  the  air.  If  the  force  caused 
by  the  deceleration  of  the  projectile  produces  stresses  beyond  the  elastic  limit  in  the  metal, 
the  projectile  may  be  deformed  or  ruptured;  this  energy  is  not  recoverable.  Similarly,  for 
the  soil,  energy  may  be  used  to  comminute,  deform  or  shear  sufficient  soil  at  the  surface 
to  form  a  crater.  This  energy  is  not  recoverable  to  rebound  the  projectile. 
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5.  Generation  of  at  least  four  kinds  of  soil  waves.  At  the  instant  of  impact,  four  types  of  soil 

waves  will  be  radiated  out  from  the  point  of  impact.  The  body  waves  which  radiate  spherically 
are  the  compression  wave  and  the  shear  wave.  The  compression  wave  oscillates  in  the 
direction  of  propagation.  The  other  two  waves  are  surface  waves  and  are  called  the  Love  wave 
and  the  Rayleigh  wave.  (Two  other  waves  called  the  coupled  wave  and  the  hydrodynamic  wave 
may  also  exist,  but  are  thought  to  be  of  minor  importance. )  The  Rayleigh  wave  causes  the  soil 
surface  to  ripple  much  like  a  wave  on  water  except  that  the  movement  of  the  particles  is  el¬ 
liptical  with  the  long  axis  perpendicular  to  the  surface  of  the  earth,  and  in  the  direction  up¬ 
ward  and  toward  the  point  of  impact.  The  Love  wave  is  a  surface  shear  wave.  The  particles 
move  in  a  path  parallel  to  the  surface  of  the  earth  and  also  at  right  angles  to  the  direction  of 
travel. 

The  velocity  of  propagation  of  these  waves  varies  approximately  as  the  square  root  of  the  density, 
and  the  period  increases  with  the  distance  from  the  origin.  The  frequency  is  believed  to  be 
approximately  proportional  to  the  velocity  cubed  divided  by  the  distance  away  from  the  point 
of  impact.  The  velocity  of  propagation  for  loose  dry  soils  is  thought  to  be  about  600  to  3,  000 
ft/sec;  for  compact,  saturated  sands  and  gravel  3,  500  to  5,  500  ft/sec;  for  clay  5,  000  to  6,  000 
ft /sec. 

Studies  perhaps  can  be  made  of  the  energy  dissipated  by  these  soil  waves;  at  least  it  should  be 
possible  from  existing  data  to  make  some  reasonable  estimate  of  the  loss. 

C.  Penetrating  Projectile 

If  the  projectile  is  rigid  enough  and  the  soil  shear  strength  sufficiently  low,  the  projectile  will  not  be 
destroyed  nor  will  it  rebound;  instead  it  will  continue  to  penetrate  the  soil  by  shearing,  comminution  or 
pulverization  and  compression  of  the  soil.  Energy  would  still  be  lost  by  the  phenomena  delineated  above. 
Energy  will  also  be  lost  if  the  soil  near  the  surface  is  sheared  and  given  sufficient  velocity  to  form  a  crater 
around  the  penetration  hole.  Nonsaturated  soils  can  certainly  be  compressed,  but  saturated  soils  are  es¬ 
sentially  incompressible  during  rapid  loading  since  the  water  does  not  have  time  to  drain  out.  The  water  and 
individual  soil  particles  are  almost  incompressible.  If  the  soil  has  sufficient  cohesive  strength,  the  penetra¬ 
tion  hole  remains  open  after  impact;  thus,  if  the  soil  is  incompressible,  the  displaced  volume  of  the  soil  once 
occupying  the  hole  must  flow  around  the  projectile  toward  the  surface  and  raise  the  level  of  the  surface. 
Energy,  of  course,  is  dissipated  in  this  operation.  For  soils  such  as  dry  sand  and  gravel  with  no  cohesive 
strength,  the  penetration  hole  probably  does  not  remain  open  and  the  surface  is  raised  by:  (1)  the  throw  out,  (2) 
the  volume  of  the  projectile  itself,  and  (3)  any  dilation  of  the  soil  caused  by  the  shearing  deformation.  If 
the  soil  is  loose  and  below  critical  density,  the  dilation  is  negative.  If  the  sand  is  saturated  and  the 
density  is  greater  than  critical,  cavitation  will  have  to  occur  between  the  soil  particles.  It  is  doubtful 
that  there  can  be  any  pulverization  of  silts  and  clays  because  these  particles  arc  already  small  and  plate - 
like.  Silt  particles  range  in  thicknesses  from  0.06  mm  to  0.002  mm  and  clay  particles  have  thicknesses 
from  0.002  mm  down  to  about  7.5  Angstroms.  If  these  soils  were  not  saturated,  compression  would  result 
from  structural  collapse.  Rock,  gravel,  and  sand  would  be  ground  into  finer  particles  by  the  energy  of 
the  impacting  projectile  if  they  are  not  saturated  and  if  the  gravel  and  sands  are  above  critical  density. 

If  those  large  grained  soils  are  comminuted,  the  void  ratio  would  increase.  The  soil  shock  waves  which 
are  initiated  when  the  projectile  hits  the  surface  will  probably  collapse  the  structure  in  the  loose,  un¬ 
saturated  materials  and  decrease  the  angle  of  internal  friction  in  all  saturated  and  in  all  dense  materials. 

In  all  cases,  the  wave  fronts  will  precede  the  projectile  by  a  very  large  distance,  since  the  projectile  even 
at  impact  is  not  traveling  at  the  speed  of  sound. 

Because  of  the  diversity  of  soils  and  the  many  different  types  of  phenomena  taking  place  during  shear¬ 
ing,  the  problem  should  be  considered  by  studying  at  least  four  different  general  types: 

Type  1.  Saturated  Loose  Sand  and  Silts  and  Saturated  Normally  Consolidated  Clay 

a.  These  materials  have  a  drained  strength  that  is  higher  than  the  undrained 
strength. 

b.  When  sheared,  the  volume  of  these  soils  tends  to  decrease. 

c.  No  drainage  can  occur  because  the  event  takes  place  in  such  a  short  length 

of  time  that  the  material  is  essentially  incompressible  and  at  constant  volume. 
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d.  The  soil  structure  is  tending  to  decrease  in  volume  upon  shearing,  and  the 
soil  is  saturated;  thus  high  compressive  porewater  pressures  result.  The 
high  pore  pressure  reduces  the  total  pressure  between  particles  and  causes 
the  effective  pressure  to  be  small. 

e.  For  a  noncohesive  soil,  the  frictional  shear  strength  in  a  plane  is  equal  to  the 
effective  pressure  on  that  plane  times  the  tangent  of  the  angle  of  internal 
friction.  If  the  effective  pressure  approaches  zero,  this  material  then  should 
behave  as  a  viscous  liquid.  This  condition  is  sometimes  referred  to  as 
'liquefaction' . 

For  cohesive  normally  consolidated  clay  with  little  frictional  strength  originally, 
the  cohesion  should  be  reduced  markedly  because  of  the  high  pore  pressure 
and  it  too  should  behave  as  a  viscous  liquid.  This  reduction  in  strength  is 
probably  a  function  of  the  moisture  content  of  the  soil. 

Of  course,  for  both  soils  the  viscous  shear  resistance  will  increase  as  the 
rate  of  strain  increases.  The  general  behavior  pattern  should  resemble  a 
non- Newtonian  liquid  and  might  be  described  by  the  Bingham  model  of  a 
non-Newtonian  liquid. 

Type  2.  Saturated  Dense  Sand  and  Silt  and  Saturated  Overconsolidated  Clays 

a.  These  materials  have  an  undrained  strength  greater  than  the  drained  strength. 

For  clays  this  is  only  strictly  true  so  long  as  the  mean  of  the  principal  stress 
is  lower  than  the  preconsolidation  load. 

b.  When  sheared,  the  volume  of  these  soils  tends  to  increase. 

c.  As  in  the  Type  1  conditions,  no  drainage  has  time  to  take  place  and  the  material 
is  essentially  at  constant  volume. 

d.  The  soil's  structure  tends  to  increase  in  volume,  but  the  only  way  this  can  happen 
is  for  voids  to  form  in  the  media.  This  condition  is  referred  to  as  soil  cavitation. 

je.  This  cavitation  causes  high  tensile  stresses  in  the  water  and  increases  the 

effective  stresses  between  the  particles;  thus  the  frictional  strength  increases 
as  the  rate  of  shear  increases  for  noncohesive  materials.  There  should  not  be 
a  very  large  increase  in  the  cohesive  strength  when  purely  cohesive  over¬ 
consolidated  clay  is  sheared  dynamically. 

In  each  of  these  soils  the  behavior  should  be  somewhat  similar  to  that  of  an  elastic 
solid,  with  the  strength  a  function  of  the  rate  and  amount  of  strain. 

Type  3.  Dense,  Dry  to  Moist  Sands,  Silts,  and  Clays 

a.  It  is  very  difficult  to  evaluate  the  initial  effective  stresses  with  these  soil 
masses  because  the  material  is  not  saturated;  therefore,  discussion  of  the  strength 
of  these  materials  is  usually  based  on  total  stresses. 

b.  When  sheared,  the  volume  tends  to  increase  and  any  water  in  the  soil  is  put  into 
additional  tension. 

c.  Because  of  the  low  moisture  content  the  clays  are  near  their  shrinkage  limits  and 
have  very  high  cohesive  strengths.  Desiccated  clay  should  behave  as  a  brittle 
elastic  solid. 

d.  Shock  waves  should  have  little  effect  on  the  cohesive  strength  of  the  shrunken  clays, 
but  the  rate  of  strain  should  influence  the  cohesive  strength.  No  locking  or  corn- 
munition  would  be  expected. 


e.  Shock  waves  in  the  dense  granular  material  should  decrease  the  strength  of 
material,  breaking  down  any  water  meniscuses  between  the  grains  and  causing 
the  particles  to  vibrate.  The  friction  strength  is  directly  proportional  to  the 
confining  pressure. 

f.  The  frictional  strength  will  consist  of  resistance  to  sliding  and  slipping  and 
resistance  to  volume  change.  The  grains  are  interlocked  and  must  roll  over 
one  another  as  shear  stresses  displace  them.  If  the  confining  pressure  is  great 
enough,  the  grains  will  be  sheared  and  broken.  This  phenomenon  is  sometimes 
referred  to  as  communition  and  gives  great  strength  to  these  materials. 

g.  These  granular  materials  then  should  behave  somewhat  like  an  elastic  solid  until 
crushing  starts.  The  crushing  phase  should  resemble  the  behavior  of  an  ideally 
plastic  solid  until  interlocking  is  once  again  brought  into  play.  For  very  high 
loads  this  process  can  repeat  itself  several  times. 

h.  The  clays  should  behave  like  an  elastic  solid  until  slippage  occurs,  then  behave 
as  a  plastic  solid. 

Type  4.  Loose,  Dry  to  Moist  Sands  and  Silts 

a.  These  materials  all  decrease  in  volume  when  sheared  and  are  compressible. 

b.  The  tension  in  any  water  in  the  soil  is  reduced. 

c.  The  fluid  within  the  voids  is  predominately  air. 

d.  Shock  waves  should  decrease  the  frictional  strength  of  sand  and  collapse  the 
structure  of  silt  which  may  be  cemented  in  its  loose  state. 

e.  The  compressive  component  of  the  stress  caused  by  the  penetrating  projectile 
as  well  as  the  shearing  stresses  can  also  cause  structural  collapse.  It  may  be 
possible  that  all  these  phenomena  will  cause  the  volume  to  decrease  to  the  critical 
volume  for  the  particular  stress  level. 

f.  If  the  structure  collapses,  the  shear  strength  may  be  relatively  unimportant;  rather, 
the  resistance  to  compression  will  be  the  most  important  characteristic  of  these 
soils  to  resist  penetration. 

g.  After  the  soil  is  compressed  in  front  of  the  projectile,  shear  resistance  will  come 
into  play  when  the  soil  reaches  its  critical  void  ratio.  This  shear  resistance 
should  be  a  constant  for  a  given  rate  of  strain  and  given  mean  principal  stress. 

D.  Proposed  Trial  and  Error  Solution 


In  the  case  of  a  soil  which  is  not  significantly  compressible,  it  may  be  possible  to  use  a  trial  and 
error  solution  to  determine  the  depth  of  penetration  of  a  projectile.  In  order  to  see  what  the  mechanism 
of  failure  is  when  a  projectile  penetrates  soil,  a  series  of  crude  model  tests  were  performed  on  bulked 
silt  and  nearly  dry  Ottawa  sand.  A  projectile  was  held  up  about  6  inches  above  the  surface  of  the  soil  and 
driven  into  the  soil  by  striking  a  rod  attached  to  the  projectile.  The  soil  was  placed  in  a  plexiglass  box 
2  inches  thick  and  2  ft.  square.  The  dry  unit  weight  of  the  sand  was  108  lbs/eu.  ft.,  and  moisture  content 
was  3  percent.  The  dry  unit  weight  of  the  silt  was  75  lbs/cu.  ft., and  the  moisture  content  about  6  percent. 
The  projectile  was  1  inch  wide,  2  inches  long  and  2  inches  deep.  The  rod  attached  was  3  feet  long.  The 
weight  of  the  steel  block  was  1. 02  lbs.,  and  the  weight  of  the  rod  was  2.  08  lbs. 

Tests  were  also  run  on  saturated  sand  and  saturated  clay  but  these  results  arc  not  shown  in  this 
report. 

Motion  pictures  were  taken  of  each  test.  The  time  increment  for  the  bulked  silt  was  2  milliseconds 
and  the  time  increment  for  the  sand  was  2.5  milliseconds.  Figures  1  through  6  show  individual  frames 
taken  of  the  entering  projectile  in  the  sand  and  in  the  silt.  The  grid  in  the  soil  is  on  2  inch  centers. 
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The  crude,  qualitative,  two-dimensional  model  tests  show  a  crater  formed  at  the  surface  and  a 
cavity  formed  around  the  projectile  throughout  penetration  for  water,  dense  sand,  and  bulked  silt.  A  cone 
of  dense  soil  forms  in  front  of  the  projectile  in  dense  sand  and  bulked  silt  and  the  rate  of  compression  of 
the  bulked  silt  can  be  observed.  There  is  apparently  no  soil  contact  on  the  sides  of  the  projectile.  It 
should  be  noted  that  the  small  thickness  of  the  soil  model  and  the  semi-rigid  plexiglass  boundaries  will 
influence  the  results  to  some  degree. 

Depending  on  the  nature  of  the  cone  formed  at  the  front  of  the  projectile,  heat  and  ablation  of  the 
reactor  core  may  or  may  not  take  place.  If  the  cone  behaves  essentially  as  a  solid,  fixed  mass,  there 
will  be  frictional  heating  only  of  soil-on-soil;  and,  in  the  brief  instant  of  the  deceleration,  none  of  this 
generated  heat  can  reach  the  vehicle. 

On  the  other  hand,  if  the  soil  originally  forming  this  cone  slides  over  the  projectile  surface  and  is 
replaced  by  other  soil  during  the  deceleration  period,  the  possibility  of  soil-on- metal  frictional  heating 
exists  with  consequent  heating  of  the  projectile. 

The  total  heat  transfer  to  a  semi-infinite  solid  may  be  expressed. 


Q  =  2  k  A  yjojira  (t^  -  t  ) 


where 


Q  =  heat  transfer,  Btu 

k  =  thermal  conductivity,  (k  =  10  Btu/hr  ft°F  for  18-8  stainless  steel) 

A  =  heat  transfer  area,  sq  feet 
9  =  elapsed  time  since  heating  began  in  hours 

a  =  thermal  diffusivity,  (a  s  k/pcp  =  .  172  ft  /hr  for  18-8  stainless  steel) 
t^  =  surface  temperature,  °F 

T  =  initial  temperature  of  vehicle  before  entering  soil,  °F 
Cp  =  constant  pressure  specific  heat,  B/lb  in  °F 

O 

p  =  density,  lb  in/ft 


For  ablation  to  take  place,  the  projectile  surface  in  contact  with  the  sliding  soil  must  exceed  a 
temperature  of  approximately  2500°F.  If  we  treat  the  projectile  as  a  semi-infinite  solid  (an  easily 
justified  assumption),  use  the  properties  of  18-8  stainless  steel,  use  a  heat  transfer  area  ten  inches 
in  diameter,  assume  that  (ti  -  tj)  is  approximately  constant  at  2500°F  and  that  the  deceleration  period 
is  in  the  order  of  one  second,  then  the  above  equation  yields  a  value  of  Q  =  500  Btu  =  389,  000  ft.  lbs. 
This  represents  over  50  percent  of  the  original  kinetic  energy  of  the  projectile  and  would  result  in  the 
removal  of  only  six  mils  from  the  face  of  the  reactor  core.  Thus,  while  the  above  figures  represent 
only  the  roughest  of  calculations,  if  future  tests  reveal  even  small  amounts  of  ablation,  it  would  appear 
safe  to  conclude  that  a  substantial  fraction  of  the  projectile's  original  kinetic  energy  is  lost  in  the  im¬ 
mediate  area  of  its  front  face.  Because  of  the  complexities  of  this  problem,  this  subject  will  not  be 
pursued  further  in  this  evaluation. 

If  the  energy  loss  by  the  five  items  listed  on  pages  9  and  10  can  be  accounted  for,  the  remaining 
energy  E  to  be  used  for  penetration  is 


E  =  Em  -  E. 


R 
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where  E  is  the  sum  of  the  energy  lost  by: 

1.  Sudden  generation  of  heat 

2.  Generation  of  sound  waves 

3.  Rise  in  temperature  of  projectile  and  soil 

4.  Compression  of  projectile  and  soil 

5.  Generation  of  four  kinds  of  soil  waves 

and  where  E  is  the  sum  of  the  kinetic  energy  plus  the  potential  energy. 

As  previously  calculated,  E^p  =  704,  371  ft.  lbs.  +  150(d)  ft.  lbs.  where  d  is  depth  of  penetration. 
Therefore,  the  remaining  energy  to  be  used  for  penetration  is  E  =  704,  362  +  150d  -  Er  ft.  lbs. 

This  energy  E  can  be  equated  to  the  force  P,  which  resists  the  projectile  at  various  levels,  times 
and  distance  traveled 

E  =  Pd. 

The  force  P  resisting  the  projectile  is  a  variable  and  is  a  function  of  the  velocity,  depth,  angle 
of  internal  friction, cohesion,  initial  effective  principal  stresses  in  soil,  and  size  of  projectile.  It  may 
perhaps  be  a  function  of  other  parameters  also.  The  resisting  force  P  for  a  given  level  of  penetration  and 
velocity  can  also  be  equated  to  the  mass  of  the  projectile  M  times  the  deceleration  a  in  the  incremental 
distance  As. 


P  =  Ma 


which  is  also 


P  = 


W 


AV 
At  ' 


W 


V2\ 


At 


where 


AV  is  the  change  in  velocity  to  velocity  during  the  change  in  time  At, 


W  is  the  weight  of  projectile,  and 


g  is  the  gravitational  constant. 

The  sum  of  incremental  losses  in  energy  AE  caused  by  penetration  in  incremental  distance  As  can 
be  written 


EAE  =  E  =  E  P  As 


where 


E  =  (704,  362  +  150d,  -  E„)  ft.  lbs. 

rt 
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and 


d  =  L  As 


Therefore, 


E  =  704, 341  +  150d  -  E 


R 


L  P  As 


■j4V4-P^) 


W 

2g 


Z  V 


1 


This  problem  can  be  solved  by  trial  and  error  if  the  resistance  to  penetration  P  can  be  computed 
for  each  given  type  of  soil  for  all  depths  of  penetration  and  velocities,  and  if  the  sum  of  energy  losses  at 
the  surface  Er  can  be  evaluated.  If  these  values  can  be  computed  the  two  unknowns  in  the  above  equation 
are  As  and  AV. 

The  velocity  Vi  at  the  surface  is  known.  If  the  resistance  at  the  surface  can  be  computed  for 
the  known  velocity,  a  small  change  in  velocity  AV}  can  be  assumed.  Then  the  incremental  depth  of  penetra¬ 
tion  Asi  can  be  computed 


P  As 
1  1 


As 


1 


W 

2gpi 


Because  AVi  =  Vj  -  V2,  the  velocity  can  be  computed  for  the  level  Asi  below  the  surface.  The 
depth  and  velocity  are  now  known.  Therefore,  the  proper  procedure  would  be:  (1)  to  compute  a  new 
resistance  to  penetration  p2  which  is  a  function  of  depth  Asi  and  V2;  (2)  then  assume  a  new  change  in 
velocity  AV2;  and  (3)  compute  AS2. 

This  procedure  could  be  repeated  until  the  change  in  velocity  AV  is  zero  or  the  final  velocity  V 
is  zero. 

In  order  to  check  the  assumptions  of  AV,  the  total  energy  lost  by  penetration  must  be  found  equal 
to  the  summation  of  the  penetration  resistance  P  times  the  incremental  distance  traveled. 


E  =  Z  P  As 


where 


Z  As  =  d 


If  the  term  of  the  resistance  to  penetration  P  times  the  incremental  depth  traveled  As  is  not  equal 
to  the  total  energy  which  can  be  lost  to  penetration,  the  procedure  would  have  to  be  repeated  with  smaller 
values  of  AV  until  the  error  is  satisfactorily  reduced. 

From  these  computations,  theoretical  time  vs  depth  of  penetration,  time  vs  velocity,  and  time  vs 
deceleration  relationships  could  be  developed  which  could  be  checked  by  model  and  full  scale  tests  for 
various  soils. 

This  whole  procedure  is  dependent  upon  computing  the  resistance  to  penetration  P  for  each  depth  and 
velocity  for  a  given  soil. 
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The  crude  model  tests  indicate  that  the  mechanism  of  failure  for  the  projectile  penetrating  soil  is 
similar  to  the  mechanism  of  failure  for  a  footing  which  is  being  loaded  slowly.  The  ultimate  resistance  of 
a  footing  may  be  calculated  many  ways  by  using: 

1.  The  theory  of  elasticity 

2.  The  classical  earth  pressure  theory 

3.  The  theory  of  plasticity  and  by 

4.  The  experimental  results 

The  best  of  these  methods  is  usually  thought  to  be  a  modification  of  the  Prandtl  solution  of  the 
bearing  capacity  problem  which  uses  the  theory  of  plasticity. 

Figure  7  is  a  sketch  of  the  deformation  pattern  and  stress  zones  after  failure  for  the  "static  bearing 
capacity  problem".  According  to  Terzaghi^  the  bearing  capacity  P  of  a  circular  footing  at  depth  d  is 
given  by  the  expression 


P  =  ttR2 


1 .  3cN  +  ydN 
c  q 


+  0.  37RN7 


where 


P  is  the  vertical  ultimate  bearing  capacity  of  horizontal  circular  footings  with  perfectly  smooth  shafts 
R  is  the  radius  of  the  footing 
c  is  the  cohesive  strength  of  the  soil 
7  is  the  effective  unit  weight  of  the  soil 
d  is  the  depth  below  the  surface  of  the  soil,  and 

Nc,  Nq,  Ny  are  bearing  capacity  factors  which  are  dimensionless  and  are  functions  of  the  angle 
of  internal  friction  and  depth. 

The  values  of  Nc  ,  Nq,  Ny,  as  given  by  Terzaghi,  are  shown  in  the  following  table. 


Bearing 

capacity 

factor 

Angle  of  internal  friction  </> 

0° 

5° 

o 

o 

15° 

20° 

25° 

30° 

35° 

0 

o 

N 

5,  7 

8 

10 

14 

19 

26 

38 

60 

100 

N 

1. 0 

1.  5 

3 

5 

9 

15 

22 

39 

85 

q 

n7 

0 

0 

0.  8 

2.5 

5.5 

12 

26 

45 

130 

1 


K. 


Terzaghi,  Theoretical  Soil  Mechanics, 


John  Wiley,  N.Y. 
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Terzaghi's  bearing  capacity  factors  were  for  foundations  near  the  surface  and  at  shallow  depths. 

This  work  has  been  extended  by  Meyerhof2  to  include  foundations  placed  at  intermediate  and  great  depths. 
Generally,  all  the  bearing  capacity  factors  increase  with  depth. 

The  bearing  capacity  equation  above  represents  an  incompressible  soil  with  boundary  conditions 
described  as  earth  pressure  at  rest,  i.e.,  the  surface  of  the  homogeneous  soil  is  horizonatl  and  extends 
infinitely  in  all  directions.  Of  course  the  lateral  and  vertical  effective  pressure  increases  with  depth. 
Therefore,  this  bearing  capacity  equation  would  not  hold  for  a  projectile  penetrating  laterally  into  a  cube. 
The  shearing  strength  of  this  soil  is  given  by  the  equation  t  =  c  +  p  tan  <f>  (Mohr-Coulomb  theory  of  rupture) 
where  c  denotes  unit  or  apparent  cohesion,  <{>  denotes  the  angle  of  internal  friction  and  p  denotes  the  effective 
normal  pressure  on  the  shear  plane. 

Resistance  to  penetration  can  be  computed  if  the  cohesive  strength  of  the  soil  and  the  bearing  capacity 
factors  are  known.  For  a  given  projectile,  the  bearing  capacity  factors  are  functions  only  of  depth  and 
angle  of  internal  friction.  Cohesive  strength  and  angle  of  internal  friction  are  functions  of  rate  of  strain; 
therefore,  it  is  necessary  to  make  some  assumptions  about  the  rate  of  strain  on  the  plane  of  failure  when 
the  pi  ojectile  penetrates  the  soil.  It  seems  reasonable  that  the  velocity  of  strain  or  strain  rate  along  the 
failure  planes  at  a  given  depth  is  about  the  same  as  the  velocity  of  the  projectile  at  that  depth.  Therefore, 
in  order  to  compute  the  resistance  to  penetration,  it  is  necessary  to  know  the  relationship:  (1)  between 
the  rate  of  strain  and  the  angle  of  internal  friction;  and  (2)  between  the  rate  of  strain  and  cohesion.  Some 
work  has  been  done  in  this  area  at  M.  I.  T.  2  and  at  Harvard.  ^  However,  additional  laboratory  work  is 
required.  The  most  promising  device  for  developing  the  required  relationship  seems  to  be  the  simple  shear 
apparatus  developed  by  H.  K.  Roscoe  at  Cambridge  University.  Other  equipment  and  approaches  can  also 
be  used. 

If  the  depth  of  penetration  computed  using  the  bearing  capacity  equations  and  the  trial  and  error  method 
outlined  above  cannot  be  verified  by  actual  field  and  model  tests,  then  other  surfaces  of  failure  must  be 
assumed  and  new  bearing  capacity  factors  computed.  Also,  if  in  future  model  tests  with  other  shapes  of 
piojectiles,  other  soils,  and  higher  initial  velocities,  it  is  observed  that  the  failure  surface  deviates  from 
that  assumed  by  Meyerhof,  Terzaghi  and  others,  new  resistance  formulas  will  have  to  be  devised  to  fit  the 
mechanism  of  failure.  One  possibility  is  to  calculate  the  shear  resistance  at  plastic  failure  for  the  pro¬ 
jectile  at  various  depths.  This  might  be  done  by  assuming  various  failure  planes  and  by  using  the  minimum 
value  of  shear  resistance  for  each  depth.  Because  there  are  an  infinite  number  of  surfaces  of  possible 
failure,  this  calculation  appears  an  overwhelming  task.  But  it  may  be  possible  by  using  a  computer  and  by 
assuming  surfaces  that  seem  similar  to  the  two-dimensional  failure  planes  shown  by  model  tests.  By  this 
same  procedure  it  may  be  possible  to  compute  the  lateral  penetration  of  projectile  moving  horizontally  into 
a  cube  of  soil,  if  the  proper  boundary  conditions  are  used. 


In  a  highly  cohesive  soil  the  unit  resistance  to  penetration  should  increase  as  the  initial  cohesive 
strength  and  the  rate  of  strain  increase.  However,  the  unit  resistance  should  only  be  slightly  influenced 
by  the  diameter  of  the  projectile  and  the  depth  of  penetration. 

In  a  cohesionless  soil  the  unit  resistance  to  penetration  should  increase  as  the  depth,  size  and 
velocity  of  projectile  and  as  initial  angle  of  internal  friction  increase. 

Thus,  for  most  soils  the  resistance  of  penetration  will  increase  with  depth.  Therefore,  the  decelera¬ 
tion  will  increase  with  depth.  This  is  just  the  opposite  for  water,  because  resistance  is  a  function  of  the 
square  of  the  velocity.  If  the  velocity  change  in  the  soil  is  linear,  the  deceleration  is  constant  and  the 
resistance  to  penetration  is  constant.  As  the  moisture  content  increases,  one  might  expect  the  resistance 
to  penetration  to  be  a  function  of  velocity  and  that  the  lower  bound  would  be  water  penetration.  A  mathe¬ 
matical  model  devised  for  saturated  soils  should  degenerate  into  the  water  entry  resistance  formulas 
when  the  mass  density  and  viscosity  of  the  soil  are  considered. 


G.  G.  Meyerhof,  "The  Ultimate  Bearing  Capacity  of  Foundations,  "Geotechnic,  V.  2,  1951,  pp.  301. 

Shearing  Resistance  of  Sands  During  Rapid  Loadings,  Report  to  Waterways  Experiment  Station  (MIT,  1962) 
_ ^rmined  Strength  of  Saturated  Clayey  Silt,  Report  to  Waterways  Experiment  Station  (MIT,  March  1963). 

A.  Casagrande  S.  D.  Wilson,  "Effect  of  Rate  of  Loading  on  the  Strength  of  Clays  and  Shales  at  Constant 
Water  Content  ,  Geotechnic,  V.2,  1951,  p.  251. 
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The  velocity  patterns  for  a  steel  block  penetrating  sand  and  bulked  silt  are  shown  in  Figures  8  and  9. 
These  data  are  very  rough  but  give  some  indication  of  what  can  be  expected. 

E.  Water  Entry  Problem 

The  water  entry  problem  should  be  studied  as  the  lower  bound  of  clays  with  very  high  saturated 
moisture  contents  such  as  swamps  and  muskegs.  Also  more  information  and  experimental  data  are  avail¬ 
able  for  this  problem  and  any  formula  devised  for  projectile  penetration  into  soils  should  degenerate  into 
formulas  for  water  entry  as  the  density  and  viscosity  decreases  to  that  of  water.  This  should  happen  as 
the  moisture  content  of  the  soil  increases.  Therefore,  a  study  of  the  liquid  stage  problem  should  result 
in  information  which  might  be  useful  in  predicting  the  action  of  a  saturated  clay  soil. 

As  clay  is  added  to  water,  the  density  and  viscosity  of  the  mixture  increases.  The  increased 
density  tends  to  increase  the  Reynolds  number  and  the  increased  viscosity  tends  to  decrease  the  Reynolds 
number  ;  consequently,  the  drag  coefficient  is  probably  not  changed  significantly,  and  the  problem  can  be 
calculated  in  the  same  way  that  the  clear-water  problem  is  treated.  Even  if  these  changes  do  not  com¬ 
pensate  each  other  or  even  if  they  both  caused  the  Reynolds  number  to  change  in  the  same  direction,  the 
drag  coefficient  seems  to  be  essentially  independent  of  the  Reynolds  number  at  values  of  the  Reynolds 
number  beyond  10^.  For  a  Reynolds  number  of  10^  and  water  as  the  medium,  the  10  inch  prototype  body 
would  require  a  velocity  of  less  than  one  foot  per  second;  but  the  terminal  velocity  in  water  is  about  18 
feet  per  second,  so  that  the  drag  coefficient  may  be  considered  to  be  a  constant  during  the  entire  de¬ 
celeration  period. 

In  Navord  report  (1308)  on  water  entry  whip  and  deceleration  of  full  scale  torpedo  models  with 
ogive  and  plate  ogive  heads,  Waugh  and  Ager  determined: 

1.  The  deceleration  at  entry  was  measured  photographically  by  analyzing  the  pictures  from  water 
entry  to  tail  submergence. 

2.  The  entry  of  the  torpedos  with  a  flat  head  was  limited  to  a  velocity  of  about  300  fps  to  350  fps 
instead  of  the  500  fps  experienced  by  the  other  shapes.  From  the  results,  the  average  entry 
decelerations  over  the  12  foot  length  of  the  torpedos  were  about  180  g's. 

3.  By  extrapolating  from  350  fps  up  to  an  entry  velocity  of  550  fps,  the  average  entry  deceleration 
over  the  length  of  the  torpedo  is  about  450  g's.  From  the  calculated  deceleration  pattern  for  the 
present  problem  involving  the  flat  nosed  cylinder,  the  acceleration  pattern  varies  from 

28,  600  fps2  at  zero  submergence  to  a  value  of  about  3000  fps2  at  a  depth  of  12  feet.  The 
average  deceleration  over  the  12  foot  length  is 


a 


avg. 


(28,  6002+  3000)  =  15  8Q0fps2 


In  terms  of  gravitational  units 


a 


avg. 


15,  800 
32.2 


=  490  g's 


This  agrees  closely  with  extrapolated  value  of  450  g's  taken  from  the  experimental  curves.  Experimental 
values  given  in  the  report  "Torpedo  Studies"  by  O.  S.  R.  D.  in  1946  are  also  comparable  to  the  values  cited 
above. 


The  Waugh  report  indicated  that  the  drag  coefficient  at  entry  is  essentially  the  same  as  though  the 
water  had  flowed  around  the  cylinder,  so  at  the  surface  the  drag  force  is  the  force  exerted  on  the  cylinder, 
except  that  the  water  surface  goes  down  and  the  force  at  entry  is  probably  not  the  force  calculated  from 
the  entry  velocity.  However,  if  compression  wave  forces  are  neglected,  the  theoretical  resisting  force 
based  on  entry  velocity  would  be: 
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p 


=  c 


D 


For  the  10  inch  cylinder  the  equation  becomes 


P  =  (0.  83)  (0.  545)  V2  =  0.  439  V2 

=  (0.439)  (550)2  =  (0.439)  (30.2)  (104) 

=  (439)  (30.2)  =  133.3  x  103  lbs. 

=  133,  300  lbs. 


where 


P  =  Drag  force  on  the  cylinder  in  lbs. 

Cp  =  Drag  coefficient 
P  =  Mass  density  of  water  in  slugs 

A  =  Cross  sectional  area  of  the  cylinder  in  square  feet. 
V  =  Velocity  of  the  cylinder  in  fps. 

Check: 

If  a  =  28,  600  fps2  and  P  =  133,  300  lbs. 


then 


P  =  Ma 


M  - 


133,  300 
28,  600 


4.  66  slugs 


where  M  is  the  mass  of  the  cylinder  in  slugs  and  a  is  the  acceleration  in  feet  per  second  squared. 


—  =  M  =  4.  66  slugs 

W  =  (4.  66)  (32.2)  =  150  lbs. 

The  average  pressure  on  the  face  of  the  cylinder  is  then: 

P _  133,  300  _ _ _  ,_  2 

A  ~  P  0.  545  =  245,  lbs/ft  =  1700  psi.  based  on  the  entire  area  of  the  cylinder. 
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Because  of  cavitation  and  separation,  the  effective  area  supporting  the  drag  force  is  probably  less  than 
the  0.545  square  feet,  so  that  pressures  in  excess  of  2000  psi  are  probably  acting  on  parts  of  the  nose 
of  the  cylinder.  The  drag  coefficient  of  0.  83  is  taken  from  Navord  Report  1308  by  Waugh  and  Ager  in 
which  they  found  the  drag  coefficient  to  be  0.  83  for  entry  velocities  of  about  350  fps.  The  coefficient  for 
the  10  inch  cylinder  in  the  present  study  should  be  about  the  same  as  for  the  torpedos  in  the  Navord 
1308  report. 

In  the  book  "Torpedo  Studies"  from  the  office  of  Scientific  Research  and  Development,  1946,  pages 
50  and  51  state  that  the  impact  pressure  is  theoretically  infinite  but,  because  of  compressibility,  the 
actual  impact  pressure  is  only 


P  =  pcVe 

p  =  (2.  0)  (4800)  (550) 

=  (4800)  (1100) 

=  5,  280,  000  lbs/ft2 
,  2 

p  =  36,  700  lbs/in 


where 


P  =  mass  density  of  sea  water  in  slugs 
c  =  velocity  of  a  compression  wave  in  water  in  fps. 

V  =  Entry  ve  ocity  of  the  projectile  in  fps. 

Although  these  elas  'c  pressures  at  entry  are  theoretically  very  high,  if  they  exist  in  reality,  they  are 
of  such  short  duratic  i  that  they  are  probably  not  highly  significant  in  the  deceleration  pattern  and  in  the 
ability  to  do  damage.  Theoretically,  the  duration  of  the  above  calculated  stress  would  be  only  a  few 
microseconds.-  Howe  -er  other  experimenters  have  found  that  there  may  be  significant  compression  wave 
pressures  occurring  at  later  stages  of  submergence.  These  effects  might  require  more  consideration 
and  study  in  this  proposed  project,  although  it  seems  very  unlikely  that  they  are  significant  in  this 
problem. 

If  the  effects  of  compressibility  are  neglected,  the  velocity  and  deceleration  patterns  can  be 
established  by  writing  the  simple  equation  of  motion  and  solving  the  deceleration  in  terms  of  the  velocity. 
This  elementary  approach  is  being  used  until  more  data  are  available,  because  all  methods  involve 
certain  unproven  assumptions.  The  basic  equation  used  here  states  that  the  drag  force  minus  the  weight 
is  equal  to  the  dynamic  force  which  is  mass  times  acceleration.  Using  a  drag  coefficient  of  0.  83  and  a 
weight  of  150  pounds  instead  of  a  submerged  weight  because  the  cylinder  is  traveling  in  a  vapor  cavity, 
the  basic  equation 


P  -  W  =  Ma 


becomes 


C 


D 


150  «  Ma 


and  then  reduces  to 
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0.094V2  -  32.2  =  a. 


This  equation  gives  the  results  in  the  accompanying  curves.  (Figure  10). 

The  calculations  show  that  the  prototype  cylinder  should  decelerate  to  a  terminal  velocity  of  18.  5 
fps  at  a  distance  of  52.8  feet  below  the  surface.  If  a  curve  of  distance  versus  time  can  be  established, 
the  slope  at  each  point  on  the  curve  represents  the  velocity.  In  turn,  the  slope  of  the  velocity  vs.  time 
curve  at  each  point  represents  the  acceleration.  On  the  basis  of  the  acceleration  curve  it  appears  that 
the  maximum  forces  acting  on  the  cylinder  would  occur  at  or  near  the  surface  of  the  water.  According 
to  the  calculations,  neglecting  compressibility,  the  approximate  force  on  the  cylinder  at  entry  would  be 
approximately  130,  000  lbs. 

To  obtain  experimental  data  for  testing  the  method  of  analysis  used,  a  small  model  was  studied  in 
the  laboratory.  A  small  steel  cylinder  weighing  0.988  lbs.  was  dropped  into  a  tank  of  water.  The 
cylinder  was  1.5  inches  in  diameter  and  2.05  inches  long.  The  entry  velocity  of  the  cylinder  as  it  entered 
the  water  was  approximately  16  feet  per  second.  Pictures  were  taken  at  400  frames  per  second  and  the 
velocity  of  the  cylinder  was  calculated  at  various  depths  of  submergence  by  photographic  techniques. 

The  results  of  this  experiment  are  plotted  in  Figure  11. 

Using  the  method  explained  above,  the  theoretical  curve  was  determined  and  plotted  for  comparison 
purposes  on  Figure  11.  The  drag  coefficient  used  was  1.1.1  The  data  from  the  experiment  are  refine¬ 
ments  not  sufficiently  consistent  to  justify  additional  trial  in  the  theoretical  procedure.  However,  the  ' 
comparison  between  the  theory  and  the  experimental  results  is  fairly  close,  indicating  that  the  theoretical 
procedure  can  probably  be  refined  and  used  to  predict  the  operation  of  the  prototype  cylinder  striking  a 
water  surface. 

F.  Computing  Penetration  in  Compressible  Soils 

For  soils  that  are  compressible,  the  Prandtl  attack  will  probably  be  inadequate  because  the 
mechanism  of  failure  is  different.  Such  nonsaturated  soils  as  sand,  and  silt  in  the  bulked  or  cemented 
state  (below  critical  density)  are  compressible  at  all  rates  of  strain  and  fail  by  compression  in  front  of 
and  lateral  to  the  direction  of  penetration.  Once  the  structure  of  the  soil  is  broken  down  and  the  density 
increases  to  critical,  then  shear  resistance  will  be  brought  into  play. 

A  trial  and  error  solution  can  be  promulgated  which  is  similar  to  the  one  proposed  in  section  C  of 
this  part,  but  compressibility  of  the  soil  and  friction  on  the  dense  cone  of  soil  on  the  face  of  the  projectile 
will  furnish  the  resistance.  From  the  crude  model  tests  (Figures  3-6),  it  can  be  seen  that  friction  on  the 
sides  of  the  projectile  is  unimportant.  The  volume  of  soil  displaced  by  the  projectile  in  a  given  time 
interval  will  be  taken  up  by  densifying  the  soil,  probably  linearly,  out  to  some  lateral  distance  to  form  a 
hemisphere  in  front  of  the  projectile. 

The  laboratory  measurement  of  dynamic  compressibility  of  soils  might  yield  the  distance  to  which 
the  soil  is  influenced  by  the  projectile  and  the  resistance  to  the  movement  of  the  projectile  at  various 
velocities.  If  so,  a  trial  and  error  solution  can  be  obtained  for  the  depth  of  penetration. 


1 


Given  on  page  3-12,  Fig.  21  of  Fluid  Dynamic  Drag  by  Hoerner. 
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VI.  Proposed  Experimental  Program 


A.  Model  Tests 


1.  Similitude  Analysis  -  Model  tests  will  be  run  to  help  study  this  problem.  In  order  that  the  results 
of  the  model  tests  represent  the  penetration  of  the  prototype,  it  is  necessary  to  consider  the  re¬ 
quirements  of  model  similitude.  In  so  doing,  scale  factors  are  derived  from  which  quantitative 
results  may  be  obtained. 

The  variables  controlling  this  modeling  problem  are  assumed  to  be: 
d  -  penetration  in  feet 

D  -  characteristic  size  of  the  projectile  in  feet 
p  -  mass  density  of  projectile  in  slugs 

V  -  velocity  of  projectile  in  feet  per  second 

V  -  volume  in  cubic  feet 

y  -  unit  weight  of  soil  in  pounds  per  cubic  foot 

a  -  deceleration  of  projectile  in  soil  in  feet  per  second  squared 

c  -  cohesion  of  soil  in  pounds  per  square  foot 

C  -  compressibility  of  soil  in  pounds  per  square  foot 

0  -  angle  of  internal  friction  of  soil  in  degrees 

g  -  gravitational  acceleration  in  feet  per  second  squared 

A  dimensional  analysis  according  to  Buckingham's  tt  -  theorem  results  in  a  functional  equation 
of  the  form: 


*1  =  f(V  V 


7 r 

n-m 


), 


where 


n  =  the  number  of  variables  involved 

m  =  the  number  of  dimensions  included  in  the  n  variables 
in  which  the  n  terms  are  sets  of  dimensionless  ratios. 

A  dimensional  analysis  (see  the  appendix)  of  the  variables  controlling  this  problem  results  in: 


d^ 

D 
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The  dimensionless  numbers  for  the  two  systems  must  be  equal  for  similitude  in  order  that  the 
unknown  function,  f,  shall  have  the  same  value  for  the  model  as  for  the  prototype.  The  sub¬ 
scripts,  m  and  p,  refer  to  the  model  and  prototype  respectively. 

Similarity  between  model  and  prototype  also  requires  that  various  ratios  between  them  be  equal. 
That  is,  all  geometric  ratios  should  be  based  to  the  same  ratio.  Thus, 


Dia.  Length 

_  _ m  _  m 

'"Ratio  ~  Dia.  Length 

P  P 


Area 

m 


Area 

P 


and 


Volume 

m 


Volume 

P 


V 

m 


V 

P 


All  force  ratios 


should  also  be  equal. 


The  ratio  of  dead-weight  forces  is  obtainedfor  the  projectile: 


W  p  g  V 
_  mm  m  m 

fr  =  nr  =  ' 

p  pP  gP  VP 


V 

m 


V 

P 


since  the  gravitational  acceleration  will  be  the  same  for  both. 
The  inertial  force  ratio  is  also  obtained  for  the  projectile: 


I 


R 


M  a 
m  m 

M  a 


P  P 


p  V  a 
m  m  m 


P  V 
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a 
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Equating  the  dead-weight  force  ratio  to  the  inertial  force  ratio  yields: 


F 


R 


R 


P  V 
m  m 


V 
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p  V  a 
m  m  m 


V 

P 


a 
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This  equality  is  true  only  if  a  =  a  .  Thus,  the  deceleration  for  the  model  must  equal  that 
of  the  prototype.  m  ^ 


Consider  next  the  functional  equation  for  the  penetration,  and  examine  first  the  ir^  term, 

7fg  =  Da/V^.  For  similarity,  the  magnitude  of  ^3  must  be  the  same  for  both  the  model  and 
the  prototype.  Therefore, 


D 

m 


a 

m 


D  a 

=  P  P 

'  v2 

P 
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And  since  a  =  a  ,  the  relationship  reduces  to: 

m  p 


D 

m 

TT~ 


p 


In  ratio  form,  this  equation  becomes: 


D 


R 


(1) 


Equating  the  7r2  term. 


*2  =  D  y/pV 


2 


for  the  model  and  prototype  yields : 


D  y  D  y 
m  m  _  p  p 


2  2 
p  V  p  V 

mm  p  p 


This  statement  may  be  rearranged  to  give  in  ratio  form: 


°R  TR  =  pr  vr  ’ 


which  may  be  simplified  according  to  Equation  1  to  yield: 


TR  =  PR  ' 


From  ^  =  c/pV 


(2) 


m 


p  V  2  p  V  2 
mm  P  P 


which  in  ratio  form,  yields: 


CR  PR  VR 


(3) 


By  similar  procedure  with  ir  =  c/pV  , 


R  PR  R  is  obtained. 


(4) 


The  sixth  ir  +  cr  m  0,  implies  that 


0 


m 


(5) 
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Now  that  the  dimensionless  numbers  are  equal,  the  unknown  function  of  the  original  equation 
for  both  the  model  and  prototype  is  equal.  Therefore,  the  ratio  of  that  equation  results  in: 


d 

m 

d 

P 


D 

m 

15“ 

p 


or  in  ratio  form: 


dR 


D 


R' 


(6) 


which  states  that  the  penetration  is  of  the  same  geometric  ratio  as  the  characteristic  size  of 
the  projectile.  To  summarize,  this  analysis  states  that  if  Equations  1  through  5  are  satisfied, 
the  model  will  be  similar  in  performance  to  the  prototype. 

2.  Description  of  Model  Tests  in  Soil  -  The  tests  will  be  begun  by  using  the  derived  equations  for 
similutude  to  establish  the  various  properties  of  the  projectile  and  the  soil  (or  other  media)  for 
the  model.  A  "model  of  a  model"  method  will  be  used  to  verify  the  validity  of  the  analysis  and 
also  to  determine  if  any  scale  effects  result  from  decreasing  the  sizes.  That  is,  the  model 
experiments  will  use  two  projectiles  (say  3/4"  and  1-1/2"  characteristic  dimension).  Thus, 
the  smaller  one  may  be  thought  of  as  a  model  of  the  larger  one.  This  is  necessary  because  of 
the  limited  amount  of  data  and  also  the  difficulty  of  collecting  that  data  for  the  full-size  prototype. 

The  models  will  be  two-dimensional,  rather  than  three-dimensional  as  in  the  prototype,  so  be¬ 
havior  of  the  soil  during  penetration  can  be  observed  with  the  Fastax  cameras.  The  projectiles 
for  the  models  will  be  of  parallelepiped  form  rather  than  the  cylindrical  form  as  in  the  prototype. 
The  width  of  the  model  projectile  will  be  compared  to  the  diameter  of  the  prototype  as  the 
characteristic  dimension. 

The  model  projectiles  will  be  shot  vertically  into  several  types  of  media  in  various  states  to 
simulate  the  prototype  entering  different  types  of  soils.  The  soil  sample  will  be  about  2  inches 
thick  and  have  a  horizontal  dimension  of  about  3  feet  and  a  variable  depth  depending  on  the  ma¬ 
terial  and  the  expected  penetration.  The  exact  size  of  the  sample  will  have  to  be  large  enough 
so  that  the  compression  shock  wave  does  not  reach  the  boundary  before  the  projectile  comes  to 
rest.  These  phenomena  will  have  to  be  observed  first  by  experiments  before  setting  the  final 
dimensions  of  the  sample.  The  container  of  the  sample  will  have  a  heavy  plexiglass  front  so 
that  pictures  can  be  made  during  the  penetration.  A  loose  grid  of  carbon  particles  will  be 
embedded  on  the  face  of  the  sample  so  that  the  movement  of  the  sample  particles  may  be  observed. 

Similar  tests  may  also  be  performed  by  shooting  the  model  projectiles  horizontally  into  a  soil 
sample.  Only  materials  with  cohesive  strength  or  artificially  imposed  cohesion  can  be  used  in 
these  tests.  From  these  simple  tests  it  will  be  seen  whether  or  not  the  present  horizontal  rocket 
tests  in  any  way  compare  to  what  will  happen  when  the  projectile  impacts  vertically  on  the  earth's 
surface. 

The  two-dimensional  tests  will  allow  several  important  observations  to  be  made: 

(1)  Mode  of  failure. 

(2)  Distortion  of  soil  below  plane  of  failure. 

(3)  Rate  of  progress  of  planes  of  failure. 

(4)  Velocity  of  particles  forming  crater. 

(5)  Depth  and  size  of  crater. 

(6)  Rate  of  deceleration  of  projectile. 
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(7)  Amount  of  elastic  rebound  of  unsheared  soil. 

(8)  Effect  of  pore  water. 

The  instrumentation  and  data  collection  will  be  very  difficult  for  such  items  as  3,  4,  6,  and  7 
but  the  net  effect  should  be  to  give  a  better  understanding  of  the  problem  and,  thus,  make  it 
easier  to  formulate  an  analytical  approach  to  it.  Three-dimensional  model  tests  will  also  have 
to  be  performed  to  determine  the  factors  involved  in  using  two-dimensional  rather  than  three- 
dimensional  tests. 

At  least  five  different  types  of  material  will  be  used  in  the  model  tests.  These  are: 

(a)  Plaster  of  Paris  to  simulate  the  non-fissured  rock. 

(b)  A  very  fine  rock  flour  to  represent  the  noncohesive  materials. 

(c)  A  highly  plastic  clay  such  as  montmorillonite  to  represent  the  cohesive  materials 
with  the  pore  fluid  either  water  or  oil. 

(d)  Some  transparent  double  refracting  material  such  as  gelatine  for  photoelastic 
studies  with  which  it  may  be  possible  to  actually  compute  the  time -pressure 
relationship  for  the  impacting  projectile. 

(e)  Water. 

In  the  case  of  (b)  and  (c),  the  density  and  moisture  content  can  be  varied  over  a  wide  range  so 
that  the  effect  of  pore  water  pressure  on  possible  liquefaction  and/or  cavitation  and  effect  of 
changes  in  shear  strength  can  be  observed. 

3.  Description  of  Model  Tests  in  Water  -  The  problem  of  the  projectile  entering  a  body  of  water 
must  also  be  considered.  In  order  to  predict  the  deceleration  pattern,  model  tests  will  be 
conducted.  The  scale  factors  for  similitude  must  be  found. 

As  the  projectile  enters  the  water,  a  cavity  is  formed.  The  size  and  shape  of  this  cavity  will 
affect  the  path  of  the  projectile.  The  Froude  number,  the  ratio  of  inertia  forces  to  gravity 
forces,  is  the  most  important  parameter  at  this  stage.  Previous  tests  by  other  investigators 
have  shown  that,  with  the  model  open  to  the  atmosphere  and  scaling  factors  derived  according 
to  the  Froude  number,  good  results  are  achieved  with  respect  to  the  cavity  and  penetration. 

The  Froude  number  criterion  requires  that: 

V  V 

m  _  p 

y°m  gm  y°p  gp 


Since  g 

m 


g  ,  the  velocity  ratio  must  be: 


VR  ^ °R  ' 


(1) 


To  find  the  time  relationship  between  the  model  and  the  prototype,  the  velocity  may  be  written 
in  terms  of  its  basic  units  and  combined  with  Eq.  (1). 
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Simplifying,  it  is  seen  that 


T 


R 


(2) 


The  geometric  ratio  between  model  and  prototype  will  remain  the  same  as  previously  discussed. 

The  model  tests  conducted  will  be  similar  to  those  conducted  with  soil  samples.  Two  different 
sized  model  projectiles  will  be  used  so  that  the  "model  of  a  model"  method  may  be  used.  From 
the  data  about  the  small  model,  predictions  of  the  velocity  and  acceleration  patterns  will  be  made 
for  the  larger  model.  The  results  from  the  larger  model  will  then  be  used  to  check  these  pre¬ 
dictions. 

Using  high  speed  photography  and  accelerometers,  the  relationship  between  depth  and  time  can 
be  established.  From  this  curve  and  from  the  photographs,  the  velocity  patterns  and  the  ac¬ 
celeration  patterns  may  be  determined.  The  areas  of  cavitation  on  the  face  of  the  projectile 
can  also  be  studied  and  determined.  From  the  accelerations,  methods  will  be  worked  out  to 
determine  the  maximum  forces  to  be  expected. 

From  the  results  of  the  two  models,  effective  means  of  predicting  the  velocity  and  acceleration 
patterns  can  be  studied.  These  studies  can  then  be  used  to  make  predictions  for  the  actual  proto¬ 
type. 

B .  Laborat ory  Tests 

Tests  on  various  types  of  earth  in  different  states  should  be  performed  to  describe  the  necessary 
empirical  constants  for  use  in  the  analytical  solution.  These  tests  should  include: 

1.  Dynamic  ultimate  shear  strength. 

2.  Dynamic  modulus  of  deformation,  that  is,  the  relation  between  shear  stress  and  volumetric 
strain,  in  the  nonlinear  elastic  range,  the  plastic  range,  and  in  the  viscous  range. 

3.  Viscosity. 

4.  Dynamic  compressibility. 

5.  Critical  void  ratio. 

In  each  case  the  following  items  must  be  considered: 

a.  Time  rate  of  strain. 

b.  Confining  pressure  or  effect  of  depth. 

c.  Pore  water  effect  on  liquefaction  and  cavitation. 

d.  Type  of  soil  with  various  void  ratios,  structure,  and  moisture  content. 

It  is  possible  that  all  of  these  effects  might  be  studied  by  using  a  modification  of  the  simple  shear 
device  invented  by  K.  H.  Roscoe  at  Cambridge  University.  A  cubical  sample  of  soil  is  confined  with 
various  pressures  and  a  deviator  stress  applied  so  that  a  state  of  uniform  shear  exists.  The  change  in 
void  ratio  and  pore  pressure  can  be  measured  as  well  as  the  strain  rate  and  load  required  to  produce  the 
strains.  The  shear  stress  can  be  applied  in  three  different  ways:  (1)  slowly,  (2)  rapidly  with  the  use  of 
explosive  charge,  or  (3)  with  a  high  speed  vibrator. 

C.  Full  Scale  Vertical  Tests 


After  the  model  tests  have  been  completed  and  some  analysis  finished,  full  scale  vertical  tests  should 
be  instigated.  Since  many  horizontal  tests  have  been  performed,  some  full  scale  vertical  tests  should  be 
performed  so  that  comparison  can  be  made  between  the  model  tests  and  full  scale  horizontal  and  vertical 
tests.  If  possible,  at  least  one  of  the  vertical  tests  should  use  a  target  constructed  with  the  same  material 
and  density  used  in  the  previous  full  scale  horizontal  tests. 
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The  targets  might  be  constructed  in  an  excavation  adjacent  to  the  tower  that  will  be  used  to  support  the 
leads  which  the  projectile  will  follow  to  the  point  of  impact.  This  excavation  should  be  about  10  feet  square 
at  the  surface  and  decrease  in  horizontal  size  to  about  4  feet  square  at  depth  of  6  feet,  so  that  the  target  ma¬ 
terial  can  be  compacted  against  the  walls  of  the  excavation  to  provide  boundary  conditions  similar  to  ones 
that  would  exist  if  the  target  were  infinite  in  size.  Precise  controls  should  be  exercised  during  construction 
of  these  targets  and  all  index  properties  determined  for  soils  used  in  the  targets.  The  layers  should  be  com¬ 
pacted  in  about  6 -inch  lifts,  so  that  the  density  and  moisture  content  is  uniform.  Evaporation  should  be  con¬ 
trolled  and  each  layer  scarified  before  the  next  lift  is  applied  so  that  no  planes  of  weakness  will  exist  in  the 
mass.  The  target  variables  will  be  size,  boundary  conditions  and  state  of  the  media.  If  the  materials  are 
compacted  into  an  excavation,  the  boundary  conditions  should  be  approximately  the  same  as  the  ones  that 
would  exist  where  the  mass  is  uniform  over  an  area  and  depth  much  larger  than  the  target.  The  variables 
which  should  be  considered  for  the  state  of  media  are: 

1.  Type,  size,  and  shape  of  mineral  particles. 

2.  Structure  of  soil. 

3.  Density  or  void  ratio. 

4.  Moisture  content. 

To  get  a  fairly  general  picture  of  the  behavior  of  the  earth's  surface  when  hit  with  a  test  item,  the 
following  type  targets  should  be  used: 

*1.  Concrete  to  simulate  hard  unfractured  rock. 

2.  Dry  loosely  placed  coarse  graded  sand-gravel  mixture.  Maximum  size  of  aggregate  1-1/2". 

3.  Dense  moist  coarse  graded  sand-gravel  mixture. 

*4.  Dense  dry  well  graded  sand. 

5.  Bulked  moist  well  graded  sand. 

*6.  Saturated  dense  sand. 

7.  Low  density  dry  cemented  silt. 

8.  Saturated  silt. 

*9.  Desiccated  medium  plastic  clay. 

*10.  Saturated  bentonitic  clay,  normally  consolidated  if  possible. 

*11.  Water. 

If  it  is  not  possible  to  perform  all  of  these  tests,  at  least  one  test  should  be  run  using  the  target  ma¬ 
terials  that  are  starred  above. 

The  variables  for  the  projectile  will  be  the  mass,  size,  velocity,  shape,  and  rigidity.  The  shape, 
size,  and  mass  should  be  limited  to  that  of  the  SNAP  reactor.  Because  of  instrumentation  difficulties,  the 
projectile  might  be  constructed  to  be  very  rigid.  The  velocity  at  impact  should  be  450,  550,  and  650  ft/sec. 
The  instrumentation  for  the  projectile  should  be  capable  of  measuring: 

1.  Rate  of  deceleration. 

2.  Pressure  distribution  over  the  face  of  projectile  versus  time,  if  possible. 

A  series  of  pressure  cells,  strain  measuring  devices,  and  accelerometers  should  be  used  to  measure 
the  effect  of  the  impacting  projectile  on  the  target.  Use  might  be  made  of  devices  now  under  development  at 
the  University  of  New  Mexico.  If  several  Fastax  cameras  are  trained  on  the  point  of  impact,  the  velocity 
of  the  expelled  material  can  be  measured.  After  the  projectile  is  at  rest,  measurements  should  be  made  so 
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that  the  penetration,  heave  of  surface,  size  of  crater,  amount  of  material  expelled,  and  amount  of  material 
that  falls  back  into  the  crater  can  be  described. 

D.  Full  Scale  Horizontal  Tests 


Even  though  many  horizontal  tests  have  been  performed,  only  one  type  of  material  has  been  used  and 
little  is  known  about  its  state.  If  the  boundary  condition  problem  can  be  solved,  a  prediction  of  behavior  of 
the  horizontal  target  and  impacting  projectile  might  be  made.  These  predictions  should  be  tested  using 
several  types  of  soil  in  different  states.  Tests  of  this  nature  would  be  useful  also  in  designing  a  target  that 
will  simulate  a  given  point  on  the  surface  of  the  earth.  By  comparing  these  horizontal  tests  with  the  vertical 
tests  using  the  same  target  materials,  it  should  be  possible  to  design  a  system  so  that  the  boundary  conditions 
are  the  same.  The  horizontal  tests  seem  to  be  very  much  favored.  Higher  velocities  and  closer  control  of 
test  conditions  can  be  obtained  on  the  long  sled  track  than  with  air  guns  mounted  in  vertical  towers.  If  the 
horizontal  rocket  sled  tests  are  to  be  of  much  use,  the  horizontal  targets  will  have  to  be  designed  so  that 
their  boundary  conditions  are  nearly  the  same  as  the  vertical  targets. 

Possibly  a  long  permanent  embankment  should  be  built  behind  and  around  the  target  area  and  each 
target  constructed  against  this  embankment.  The  resistance  of  this  embankment  with  proper  vertical  sur¬ 
charges  on  top  of  the  target  might  approximate  the  vertical  and  horizontal  resistance  of  the  earth. 

The  soil  target  might  also  be  constructed  between  two  steel  bulkheads  and  pressure  applied  laterally. 
The  pressure  should  increase  from  zero  at  the  face  to  some  maximum  value  at  the  rear  of  the  target  to 
simulate  the  increasing  confining  pressure  of  the  earth.  A  bulkhead  placed  on  top  of  the  target  with  pressure 
applied  in  the  same  way  would  complete  the  boundary  conditions. 

In  some  way  the  face  of  the  target  must  be  restrained  so  that  it  is  not  strained  by  slope  instability 
because  of  the  weight  of  the  material.  Whatever  is  used  will  have  to  be  expendable  and  not  interfere  with 
the  penetration  of  the  projectile. 

Essentially  this  will  be  a  trial  and  error  procedure.  A  design  must  be  generated  and  tested  until  the 
results  are  the  same  as  the  ones  obtained  with  the  vertical  tests.  It  may  very  well  be  that  different  designs 
will  have  to  be  used  for  different  types  of  material. 

It  might  also  be  possible  to  construct  a  counterbalanced  system  in  which  the  soil  target  is  rotated  in 
the  plane  of  the  sled  track  and  the  vertical  direction,  so  that  the  centrifugal  force  is  similar  to  the  force  of 
gravity.  The  boundary  conditions  imposed  could  be  made  almost  identical  to  those  on  soil  at  rest.  Of  course 
it  would  be  difficult  to  synchronize  the  rotation  of  the  target  and  the  movement  of  the  rocket  sled. 

Once  the  designs  are  completed,  specifications  should  be  prepared  for  the  construction  of  future  targets 
which  will  simulate  behavior  of  various  portions  of  the  earth  when  different  types  of  projectiles  are  impacted 
at  different  velocities. 


VII.  Survey  of  Similar  Work  Being  Done  Elsewhere 


At  the  outset  of  the  Sandia  Earth  Target  Simulation  Study,  it  was  recognized  that  diligent  efforts  must 
be  exercised  to  insure  that  the  work  proposed  to  be  done  in  this  study  will  not  repeat  or  duplicate  work  being 
done  by  others.  At  the  same  time,  it  was  recognized  that  a  close  survey  of  other  work  being  performed  in 
the  soil  mechanics  field,  even  though  not  directly  connected  to  this  study  or  its  objectives,  might  yield  very 
useful  data,  ideas,  equipment,  and  techniques  which  could  be  profitably  used  in  Sandia's  study.  Accordingly, 
conferences  were  held  with  several  laboratories  engaged  in  soil  mechanics  research  related  to  the  problems 
under  study. 

Several  institutions  are  researching  the  basic  behavior  of  soil  under  dynamic  loadings,  but  most  of 
this  work  is  directed  toward  blast  effects  on  underground  structures,  soil  wave  propagation,  and  the  effect 
of  vibratory  and  transient  loads  on  individual  footings.  Of  particular  interest  to  this  program  are  the  stress- 
strain-time  relationships  that  have  been  studied  for  various  soils.  The  type  of  equipment  (and  the  limitations 
of  it)  used  in  these  tests  is  also  of  interest. 
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Some  of  the  institutions  and  the  areas  in  which  they  are  working  are  listed  below: 

1.  Massachusetts  Institute  of  Technology,  Cambridge,  Mass. 

a.  Dynamic  Triaxial  Strength  Tests 

b.  Seismic  velocity  in  confined  samples 

c.  Soil  structure  interaction 

2.  University  of  Illinois,  Urbana,  Illinois. 

a.  Footings  subjected  to  transient  loads 

b. .  Energy  absorption  of  granular  materials 

c.  Soil  structure  interaction 

3.  U.  S.  Corps  of  Engineers  Waterways  Experiment  Station,  Vicksburg,  Miss. 

a.  Footings  subjected  to  transient  and  vibratory  loads 

b.  Wave  propagation  in  confined  samples 

c.  Dynamic  soil  structure  interaction 

4.  U.  S.  Naval  Civil  Engineering  Laboratory,  Port  Hueneme,  Calif. 

a.  Wave  propagation  in  confined  samples 

b.  Footings  subjected  to  blast  loads 

5.  Illinois  Institute  of  Technology  Research  Institute,  Chicago,  Illinois. 

a.  Wave  propagation  in  confined  and  triaxial  samples 

b.  Footings  subjected  to  transient  loads 

c.  Dynamic  soil-structure  interaction 

d.  Dynamic  photoelasticity 

6.  University  of  New  Mexico,  Air  Force  Shock  Tube  Laboratory,  Albuquerque,  New  Mexico. 

a.  Dynamic  soil-structure  interaction  studies 

b.  Dynamic  triaxial  tests 

c.  Dynamic  modulus  of  deformation  in  the  constrained  state 

d.  Dynamic  shear  modulus 

e.  Wave  propagation  in  cylindrical  samples 

f.  Static  arching  studies 

7.  Sandia  Corporation,  Albuquerque,  New  Mexico. 

a.  Theoretical  missile  penetration  studies 

b.  Shock  waves  in  solids 


c.  Cratering  in  soils 


The  above  list  certainly  does  not  include  all  the  groups  working  in  this  area.  There  are  some  industrial 
and  consulting  engineering  firms  as  well  as  other  universities  doing  excellent  work,  but  the  general  pattern 
of  subjects  being  investigated  follows  those  listed  above. 

The  work  concerning  transient  loads  on  footings  is  somewhat  applicable  to  this  study.  All  of  this  work 
so  far  reported  deals  with  a  plate  placed  on  a  confined  soil  and  then  given  an  impulse  or  blow  by  some  sort 
of  an  actuator  or  a  blast  on  the  soil  surface  and  plate.  The  surface  deformation  and  depth  of  settlement  is 
reported.  The  soil  is  usually  dry  Ottawa  standard  sand.  This,  of  course,  is  quite  a  different  problem  from 
a  projectile  moving  at  high  velocity  and  impacting  on  the  soil,  but  some  of  the  failure  mechanisms  are 
similar. 

The  dynamic  triaxial  strength  tests  as  related  to  rate  of  strain  are  interesting  but,  like  all  triaxial 
tests,  the  strains  and  stresses  are  not  uniform  throughout  the  sample  and  it  is  very  difficult  to  decipher  what 
is  truly  happening  to  the  soil  on  the  plane  of  failure  at  the  time  of  failure. 


VIII.  Time  Schedule  for  Performance 


A  study  such  as  described  in  this  report  is  quite  difficult  to  time -schedule  accurately.  It  is  very 
desirable  to  have  such  a  time-schedule  to  predict  the  expected  duration  of  various  phases  of  the  analytical 
and  experimental  phases  of  the  entire  study.  However,  it  must  be  recognized  that  the  analysis  and  deter¬ 
mination  of  phenomena  relating  to  the  fundamental  mechanisms  of  interaction  between  nonhomogeneous 
physical  mediums  is  sometimes  difficult  and  laborious,  and  time -estimates  are  sometimes  not  met. 

In  the  light  of  the  above  statement,  a  proposed  time -schedule  for  the  performance  of  the  analytical 
and  experimental  phases  of  the  earth  target  simulation  study  has  been  made  and  is  shown  on  Figure  13. 


IX.  Conclusions 


Several  problems  have  existed  for  a  long  time: 

1.  Determining  the  depth  of  penetration  of  a  projectile  impacting  into  the  earth, 

2.  The  fundamental  failure  mechanism  involved  in  earth  impacts  of  projectiles,  and 

3.  Simulating  such  earth  impacts  with  a  horizontally  impacted  earth  target. 

Through  the  methods  presented  in  the  proposed  study,  it  may  be  possible  to: 

1.  Evaluate  the  present  horizontal  earth  target's  ability  to  simulate  the  behavior  of  the  earth 
during  impact. 

2.  Design  horizontal  earth  targets  that  will  simulate  the  various  portions  of  the  earth. 

The  proposed  research  may  yield  a  mathematical  model  which  will  describe  the  behavior  of  the  soil 
during  impact.  This  model  will  depend  on  whether  there  exists  a  unique  relationship  between  void  ratio, 
mean  principal  stresses,  shear  stress,  and  rate  of  strain  for  the  various  states  of  a  given  material.  There 
is  strong  evidence  that  such  relationship  does  exist.  These  soil  parameters  are  very  difficult  to  measure 
at  high  rates  of  strain,  but  with  proper  equipment  this  should  be  possible.  However,  it  is  also  possible  that 
the  mathematical  model  may  be  so  complicated  that  only  a  skilled  applied  mathematician  can  solve  it  with 
the  aid  of  a  computer. 
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Figure  1.  Preliminary  Laboratory  Test  of  Penetration  Into  Sand 


Figure  2.  Preliminary  Laboratory  Test  of  Penetration  Into  Sand 


Figure  3.  Preliminary  Laboratory  Test  of  Penetration  Into  Sand 
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Figure  5.  Preliminary  Laboratory  Test  of  Penetration  Into  Bulked  Silt 
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Figure  /.  Preliminary  Laboratory  Test  of  Penetration  Into  Water 


Approaches  25,  000  fps  at  t  =  0.0015  sec. 
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Figure  10.  Theoretical  Displacement  -  Velocity  -  Acceleration  of  Water  Entry  of  Snap  10A  Core  Vessel 


Figure  11  Velocity  Pattern  for  Cylinder 
Dropped  into  Tank  of  Water 
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Figure  11.  Velocity  Pattern  for  Cylinder  Dropped  into  Tank  of  Water 


after  failure 


Figure  12.  Cross-Section  Illustrating  Prandtl  Plastic- Equilibrium 
Theory  for  Footing  Bearing  Capacity 
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APPENDIX  A 


Variables  controlling  penetration: 


Symbol 

Dimensions 

Description 

d 

L 

penetration 

D 

L 

size  of  projectile 

P 

ft2l"4 

density  of  projectile 

V 

lt"  1 

velocity  of  projectile 

a 

lt"2 

deceleration  of  projectile  in  soil 

7 

fl'3 

unit  weight  of  soil 

c 

fl-2 

cohesion  of  soil 

C 

fl"2 

compressibility  of  soil 

0 

none 

angle  of  internal  friction 

The  assumption  is  made  that  the  penetration  of  a  projectile  impacting  the  earth's  surface  is  dependent 
on  the  variables  listed  above  so  that 

d  =  f(D,  p,  V,  a,  7,  c,  C,  0). 

This  equation  states  that  the  dependent  variable  d  is  a  function  of  the  eight  independent  variables  as  shown. 
Thus,  the  problem  is  one  dealing  with  nine  variables. 

In  order  to  reduce  the  number  of  variables  involved,  the  Buckingham  Pi  Theorem  is  applied.  This 
procedure  reduces  the  problem  from  one  dealing  with  nine  dimensional  variables  to  one  dealing  with  six 
dimensionless  variables.  Since  the  nine  dimensional  variables  involve  three  dimensions,  the  dimensional 
variables  can  be  reduced  from  nine  to  9-3  =  6  dimensionless  variables  such  that 

*1  af<w  wV- 

This  equation  means  that  the  jt  term  involving  the  dependent  variable  is  some  function  of  the  five  7r 
terms  which  involve  independent  variables.  An  analysis  of  this  type  often  helps  to  organize  a  research 
project  more  systematically  and  also  allows  an  experimenter  to  deal  with  dimensionless  variables  so  that 
the  number  of  variables  involved  is  smaller  than  for  the  dimensional  case. 

The  Pi-Theorem  operates  as  illustrated  in  the  following  development: 

where 

L  =  Dimension  of  length 
F  =  Dimension  of  force 
T  =  Dimension  of  time 
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The  functional  equation  is: 


d  =  f(D,  p,  V,  a,  y,  c,  C,  0). 

Choose  as  the  repeating  variables  D,  p,  V  so  that  all  three  dimensions  are  included  in  these  variables. 
Then  the  repeating  variables  are  combined  with  each  of  the  remaining  variables  in  the  following  way: 

X1  Y1  Z1  1 
n  =  D  p  V  d 

where  1  is  an  arbitrary  number  selected  as  the  exponent  of  d  or  in  dimensional  form: 

L  1(FT2l'4)  1  (LT-1)  1L1  =  L°F°T° 

Since  each  ir  term  is  to  be  dimensionless,  the  sum  of  the  exponents  of  each  dimension  must  equal  0. 

L:  X  -  4Y1  +  Z1  +  1  =  0 

F:  Y  =  0 

T:  2Y  -  Z  =  0  .  '  .  Zl  =  0 

from  the  equation  for  L 

X!  =  -l 


and  substituting  the  exponents  into  the  initial  equation 


d 

*1*5 


X2  Y2Z2  1 
7T2  =  D  p  V  7 


(a) 


L  2(FT2L'4)  2(LT_1)  2(FL  3)  =  L°F°T0 


L:  Z2  '  4Y2  +  Z2  -  3  =  ° 


F:  Y2  +  1  =  0  .  Y2  =  -1 

T:  2Y2  -Z2  =0  Z2  =  -2 

X2  =  1  from  the  equation  for  L 

and 


IT 


2 


(b) 
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X3  Y3  Z3  1 
ir  ^  =  D  p  V  a 


L  3(FT2l'4)  3(LT_1)  3(LT'2)  =  L°F°T° 


L:  X3-4Y3+Z3  +  1=0 


F:  Y3  =  0 


T:  2Y3  -  Z3  -  2  =  0  .  '  .  Zg  =  -2 


Xg  =  1  from  the  equation  for  L 


Da 

^3  =  v2 


X.  Y,  z„  , 
„  4  4  4  1 

ir  =  D  p  V  c 
4 


L  4(FT2L"4)  4(LT  *)  4<FL'2)  =  L°F0T° 


L:  X.  -  4Y,  +  Z,  -2  =  0 
4  4  4 


F:  Y.  +  1  =  0  .  '  .  Y  =  -1 

4  4 


T:  2Y  -  Z  =  0  .  ’  .  Z  =  -2 

4  4  4 


X^  =  0  from  the  equation  for  L 


c 

*4  =  T 
pV 


c 

pV2 


and  since  0  is  dimensionless  already,  p,  V,  and  D  will  disappear  and 


Similarly, 


*6  =  * 
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The  functional  equation  is  now: 


*1  “f(W  wV 


or 


d_ 

D 


0 
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